Intermittent dissipation of a passive scalar in turbulence 
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Probability density function (PDF) of passive scalar dissipation 'P(e) is found analytically in the 
limit of large Peclet and Prandtl numbers (Batchelor-Kraichnan regime) in two dimensions. The 
tail of PDF at e ^ (e) is shown to be stretched exponent In'P(e) oc e^^'"'' . At e <C (e), P oc l/\/e. 
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Probability distribution of the gradients of turbulent 
fields is probably the most remarkable manifestation of 
the intermittency of developed turbulence and related 
strong non-Gaussianity. A typical plot of the logarithm 
of gradient's PDF (which would be parabolic for Gaus- 
sian statistics) is concave rather than convex, with a 
strong central peak and slowly decaying tails. This is 
natural for an intermittent field since rare strong fluc- 
tuations are responsible for the tails while large quiet 
regions are related to the central peak. In particular, 
such PDFs were observed for the dissipation field (square 
gradient) of passive scalar advected by incompressible 
turbulence which is the subject of the present paper. 
We consider scalar advcction within the framework of 
the Kraichnan model assuming velocity field to be delta- 
correlated in time ^ . Most of the rigorous results on tur- 
bulent mixing has been obtained so far with the help of 
that model which is likely to play in turbulence the role 
Ising model played in critical phenomena. High-order 
moments of the scalar were treated hitherto by the per- 
turbation theory around Gaussian limits. Clearly, the 
kind of strongly non-Gaussian PDF observed for gradi- 
ents cannot be treated by any perturbation theory that 
starts from a Gaussian statistics as zero approximation. 

Since we consider developed turbulence with large 
Peclet number Pe (measuring relative strength of advec- 
tion with respect to diffusion at the pumping scale), it 
is tempting to use Pe^^ as a small parameter. Yet any 
attempt to treat diffusion term perturbatively is doomed 
to fail because the PDF of the dissipation is a nonpertur- 
bative object with respect to the inverse Peclet number: 
it is zero at e 7^ and zero diffusivity yet have nonzero 
limits as diffusivity goes to -f 0. Indeed, the main contri- 
bution into dissipation is given by the scales around the 
diffusion scale where advection and diffusion are compa- 
rable. Still, the presence of a small parameter calls for 
finding a proper way to simplify the description. Follow- 
ing 0-^, we show here that the dynamical formalism of 
Lagrangian trajectories is a proper language to describe 
the probability distribution of the dissipation field. Our 
goal is to express unknown (statistics of dissipation) via 
known (statistics of pumping). Since Peclet number is 



the ratio between pumping and diffusion scales, then any 
piece of scalar has a long way to go between birth and 
death, our goal is to describe how statistics is modified 
along the way. Dynamical formalism explicitly reveals 
the presence of two different time scales, a short one re- 
lated to stretching and a long one related to diffusion 
(which eventually restricts the process of stretching) ; this 
time separation has been exploited first in solving simi- 
lar problem for one-dimensional compressible velocity . 
Time scale of stretching fluctuations is of order of inverse 
Lyapunov exponent while the whole time of stretching is 
InPe times larger. Clear time separation in the dynami- 
cal history is guaranteed for (but probably not restricted 
by) trajectories that correspond to the value of e larger 
than that produced by the pumping, (e) j Pe^ , yet smaller 
than (e) InPe (whether it is possible to extend the time- 
separation procedure for a wider interval of e will be the 
subject of future analysis [^). At Pe — > cx), we are able 
to calculate PDF V(t) rigorously, exploiting time sepa- 
ration and executing explicitly separate averaging over 
slow and fast degrees of freedom. 

There are three main steps in deriving V(t) : 

1. Presenting the PDF as an average of a functional of 
the time-ordered exponent of the strain matrix; 

2. Reparametrization of the problem into nonlocal aux- 
iliary quantum mechanics 

3. Implementing time separation which makes reduc- 
tion to two coupled yet local quantum mechanics: one- 
dimensional one describing the long evolution and two- 
dimensional one describing the fast fluctuations. 

Let us define the problem and make the first step of 
the derivation. Advection of a passive scalar 9{t, r) by an 
incompressible flow v(t, r) is governed by the equation 



{dt + VaV a - kA)9 = 4> , VqWq=0, 



(1) 



where ip{t, r) is the external source and k is the diffu- 
sivity. Both v(i, r) and 0(i, r) are Gaussian independent 
random functions, ^-correlated in time |l- 



(2) 



(0(ti,ri)0(t2,r2)) = S{ti - t2)x{ri2) 



Here x(''i2) as a function of ri2 = |ri — decays on the 
scale L and x(0) is the production rate of 9^. We con- 
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sider Prandtl number Pr = vjn, (viscosity-to-diffusivity 
ratio) to be large so that L lies in the viscous interval 
where the isotropic velocity statistics is described by the 
pair correlation function {va{t\^V\)vp(ti^V2)) given by 



<5(ii-t2) 



(3) 



for the scales less than the velocity infrared cut-off L„, 
which is supposed to be the largest scale of the problem, 
represents two first terms of the expansion of the ve- 
locity correlation function in so that D ^ Vq/L^. 
We presume also the inequality Pe^ = dDL^/2K 3> 1 
which guarantees that the mean diffusion scale = 
2y/ k/D is much less than the pumping scale L. It follows 
from (^) that the correlation functions of the strain field 
(Tq/3 = V/3Va are r-independent. That property means 
that (7c(/3 can be treated as a random function of time t 
only. To exploit that, it is convenient to pass into the co- 
moving reference frame that is to the frame moviiig with 
the velocity of a Lagrangian particle of the fluid : 



dt9 + r,a,i{t)die - KAe ^ (j){t,r), 
Making spatial Fourier transform and introducing 



W(t) = Texp / (T(r)dT, 
Jo 

one can write the solution of (0) as follows 



(4) 



(5) 



X exp 







W-\t)W-'-'' (r) 



drku 



fiiy 



(6) 



Averaging over a simultaneous product of 2n-th repli- 
cas of the inverse Fourier transform of kS^, we get the 
n-th moment of the dissipation field e = niyO)^. The 
PDF restored from all the moments is 
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dt / dqxq 



qW{t)m 
27rPe 



qA(t)q' 



exp 



Pe2 



(7) 

(8) 
(9) 



where q = kL and an extra integration over the auxil- 
iary vector field m takes care of combinatorics and sum- 
mation over vector indices. Notice that the direction of 
time in (^;^ is opposite to that in (0). The averaging 

a b + c 
b — c —a 



A{t) = DW{t) / dt'W-\t')W~^^' (t')W^ it), 



over the traceless a = 



is according to 



Va{t) exp[- dt'ia"^ + b"^ + c^/2)/2D]. The whole ex- 
pression (^) is invariant with respect to the global (time 
independent) rotation m Rm, a R&R^, provided 
the pumping is statistically isotropic. That allows to get 
rid of the angular integration in dm, counting all the 
dynamical angles from the direction of m. 

Let us do the step 2 now. The general tool to aver- 
age a function of W, particularly (^), is Kolokolov trans- 
formation replacing T exp by a regular function of new 
fields 0]. We represent W = R^Vif, where R^ is the 
matrix of rotation by the angle dt'd{t')/2, V is di- 
agonal matrix with the elements exp [± /p dt'ait')] and 

Q "^^^^ ^, ip{0) = to provide for W{0) = i. 

Our transform {a,b,c} {a,'d,ip} is a hybrid version 
of those used in [||J^ and in Appendix B of The re- 
lation between the old and new variables is obtained by 
differentiating R^V(p and comparing with dW /dt = aW: 





t 




a + ib — a exp 


i Jdt'd 


-f i-exp 




. 





dt' (2(7 + iz9) 



2c = — -|- ipexp 



2 I dt'a 





(10) 



The Jacobian of the transformation J ~ JuiJi has a 

rt 



standard product Jul 
trivial Jacobian J; 



exp 



nLo exp 

J^dta{t) 



2Qdt'a 



and non- 



associated with the 



positive Lyapunov exponent, describing the exponential 
stretching of trajectories. The rotation matrix i?^, which 
is the same for all the dynamical matrix processes, can be 
removed by collective transformation of all the external 
vectors in the problem i.e. one may explicitly integrate 
over since we average i9-independent objects in (|^). 
Averaging some function of W{T) is then reduced to av- 
eraging the same function of V{T)(p with respect to the 
measure T>a{t)'Dip{t) Y\i^Qexp[2 dt'cr] exp[— 5], where 



S = 



1 
2D 



dtla^- 2Da 



exp 



4 J dta 
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(11) 



Finally, we apply the time-separation procedure to 
capture the term in (^ dominant at large Pe. Two very 
different time scales are clearly seen in (|7[-|) at least for 
e ^ (e)/Pe^ when the integral over time in (||) has to 
be large. Large value of the integral is achieved on such 
realizations of a where the integrand grows exponentially 
due to yV(t) until the last exponential factor restricts the 
grows at the time of order D^^ InPe. Well before, on a 
time scale of the inverse stretching rate D~^, the expo- 
nentially growing yV(t') makes the integral over dt' in (|) 
saturated. 
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A dynamical picture in r-space would be as follows: 
For the Lagrangian trajectories giving the main contribu- 
tion into (e"), fluid particles start very close. Diffusion re- 
mains dominant until the particles separate by a distance 
comparable to the diffusion scale r^. This phase of the 
dynamics takes place on times of order D^^, notice that 
this time is diffusion-independent since is a scale where 
diffusion and stretching are comparable. Once the dis- 
tance between particles outgrows r^, random multiplica- 
tive stretching due to velocity becomes dominant. Due 
to a multiplicative nature of the dynamics, the time to go 
from rd to L is proportional to D^^ InPe. Let us now in- 
troduce a separation time to satisfying 1 ^ Dto ^ ln[Pe], 
that time will disappear from the final answer. 

The temporal separation makes it possible to substi- 
tute t by to as an upper limit in integration over t' in (|^). 
Also, an important manifestation of time separation is a 
distinctively different behavior of the (p field (responsible 
for the "rotation" in the pseudospace) at time intervals 
smaller and larger than the separation time. The action 
(|ll|) shows that for t > to ^ 1/D the (p field is frozen 
- no dynamics at all, (p{t) = ip{to), which has a clear 
physical meaning since stretching proceeds along one di- 
rection. At the smallest times the v^-field dynamics is 
essential. Indeed, let us denote p{t) = a{t')dt' and 
explicitly transform A(i) according to (|lC|): 



A(t) 



=2^ 



\\D j dt' exp[2p{t') - 2p{t)] 



(12) 



Pit) - pit') 1 A 1 



For trajectories with predominantly positive a we replace 
pit) by pito) and neglect the integral from to to in Q at 
t > to- Since [p{t) - (^(t')]^ exp[2p(t')] and exp[-2p(t')] 
decrease exponentially as t' increases, then the integral 
over dt' saturates, which allows for time separation too. 
Finally, non-diagonal and lower diagonal elements in A 
are exponentially small. For time separation to be com- 
plete, the integral pit) at t ^ to has to be counted from to 
i.e. we replace p(t) —f pit) + pito), then in the dominant 
order in Peclet we get 

oo 

^-(2^/^^/ '^qX,.?e^^exp(-,?e^^-g), (13) 

to 

p = D r dt'e2''(*'){e-4p(t') + [^it) ^ pit')]''}, (14) 
Jo 

where we denoted P — exp[2p(to)]. Next, one needs 
to average exp [— sQ] over the short-time Va^Vp and 
a long-time "Day separately. The corresponding weights 
of averaging with respect to (t> and (t< are completely 
decomposed: 5 = 5< +5>. The great advantage of ( p^ ) 
is that, in the long-time averaging, both (3 and /i arc just 
external parameters, depending neither on time t nor on 



(T> . Once the average over cr> is performed, we are left 
with a function of fjL. The final result for P(e) will be 
then obtained by averaging over (t< and p. 

We start from averaging exp[— sQ] over (T> which re- 
duces to the following auxiliary quantum mechanics 



V>im,p) = 



,-(T-to)/2 



$(r-to;p) 



p=0 



(9t-^>) 



$ ^ 0. Hy 



■/.qx,.?e-exp 



2 P 



(15) 
(16) 



2 2p(t)_^ 



(27r)2 Pe^ 

Here sQ, with Q from (O], plays the role of potential, 
the quadratic in a term of (|l l|) gives the derivative part of 
the Hamiltonian, while the linear in a part of ( pT| ) gives 
the initial condition at t = for the "wave function": 
<i>(0;p) = e''. Since the "potential" part of the quantum 
mechanics vanishes at p — > cx) and $(0; p) does not, one 
obtains the asymptotic behavior of $ at t — > oo 



1 

2^ 



„t/2 



yPe$4zj), V> =<^.iiVcpiP)-^) 



^dyy^dy + C/(y))$,(2/)=0, 



Uix)^ 



i2n) 



j dkxkki exp [—kfx'^ 



(17) 
(18) 

(19) 



where the new variable y — e''(Pe/i/3)^^ has been intro- 
duced. <l>*^latj/— s-oo and should vanish at 
y — > 0. The generating function is a function of a sin- 
gle argument sm'^/p, we denote = Ci2sm? / p). The 
potential U (x) is everywhere positive, it has to turn into 
a constant at a; — > and behave as at large x. Un- 
der such a general assumption on the pumping function 
Xq, one can show that the only singularities of C(z) are 
poles on the negative semiaxis at a finite distance from 
zero (see Q for the details). The asymptotic condition on 
<&*(?;) gives the normalization C (0) = 1 and dX z ^ oo 
ln[l/C(z)] ~ ^/z. For example, in the particular case 



Ucix) - 



1, 
l/a;3 



a; < 1, 
X > I. 



(20) 



one finds 2C(z) = z [I^ (2Vi)/i iv^) + h (2 V^)/( ( V^)] ^ 
and C(z) iz^/'/An) exp [-3^^ at z ^ -|-oo. 

Now we have to average — C(2sm^//i) over cr< and 
p, what is equivalent to averaging over the random vari- 
able /i (responsible for the fluctuations of the diffusion 
scale due to strain variations): (e~'*'^)o- — (C(2sTO^//i))p. 
A convenient auxiliary object to calculate is the generat- 
ing function Vx — (exp [— A/i]), which can be found from 
yet another quantum mechanics: 

- / dp^ ip = 0, ^) , (h< + i") * = 0, (21) 



f 

J —< 



dpp"-^e~P 
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iJ< ^ -2e-^% - \dl + Ae-2P (i + ^^^^p) , 

where the static equation ( pT| ) follows at to ^ 1 from the 
respective dynamical one similarly to the way ( |l7| ) fol- 
lows from Exact solution of (|l]), ^-tt^I + Lp^e-^P = 



7'(m) = (2V) '^'exp(-l/[2A.]) 



(22) 



Note that ^"^/^ (which can be interpreted as inverse lo- 
cal diffusion scale) has exactly Gaussian PDF, that seems 
to be a consequence of strain Gaussianity. Integration of 
Ci^sra^ I [l) with V{[l) gives the final answer: 



P(6) 



1 
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dsC (is) / dxexp isx 



2 \/e 



(23) 



That formula is our main result. It expresses dissipation 
PDF in terms of the function C determined by ( |l8| ) with 
the potential (^9|) given by the pumping. For any pump- 
ing, all is necessary to get V{€) is to solve an ODE of 
the second order. In particular, one can find C(z) as a 
series in z iterating 4" in (|l^) over U. That expresses 
dissipation moments directly via the pumping 



(e"> = 2" [(2n- 1)!!] n! x 

oo Vi ( 

yf J 





dy2ytU{y2) ■ 



(24) 

dy2ny2nU(y2n)- 



y2n-l 

V2n-2 



One gets, particularly, the conservation law identity (e) — 
— 2C"(0) = x[0]/2 as a consistency check of our pro- 
cedure. The form of the PDF at e ~ (e) is pumping- 
dependent while the asymptotics at small and large e 
have universal forms. 

At small e (and infinite Pe) the two first terms have 
the form Vie) Ae"^/^ Bine so that gradient's PDF 
V{Ve) ^ A- B\S70\ In |V6'|-i tends to a constant (equi- 
librium cquipartition) , notice nonanalyticity of the first 
non-equilibrium correction. The asymptotics is sensitive 
to the fast angular dynamics. 

Since the only singularities of C{iz) are poles on the 
imaginary axis then large-e asymptotics is determined 
by the pole nearest to the origin in the upper half-plane: 
z = iT*,T* - 1. In particular, « 2.7 for (|2^). The 
dx-integration in (03) is done by a saddle-point method: 



Vie) 



-1/2 



exp 



(-3(27 



Nl/3 



/2 



(25) 



Stretched-exponential tail is natural for steady PDF of 
the gradients |^ contrary to lognormal unsteady distri- 
bution which takes place without diffusion Q . The value 
1/3 for the stretched exponent may be explained as fol- 
lows: Since the local value of VO is proportional to the 



scalar fluctuation at the diffusion scale 59 (which has 
exponential PDF tail |,|,0) times an inverse local dif- 
fusion scale [which is Gaussian according to (^2|)1 then 



which corresponds to (E 



Note that the asymptotics (|25|) is determined by the dy- 
namics of stretching (not of rotations), thus it is likely 
to take place in any dimensions. The data of numerics 
done directly for Kraichnan model are satisfactory 
fitted by 1/3 stretched exponent within the window ex- 
panding when Pe grows. The exponent 1/3 agrees also 
with the values 0.3-;- 0.36 given by numerics ^ and 0.37 
by experimental data on air turbulence p^ . Moreover, 
our exponent derived formally at Pr ^ 1 agree with the 
data of 13 1 obtained at Pr ~ 1 as well. This is sur- 
prising because 7^(e) is independent of Pe at large Pr 
while at Pr < 1 the statistics of e is markedly different 
since (e") oc Pe^" where /i„ are anomalous exponents of 
the scalar ||lj]. That may be, however, that the PDF tail 
is still given by ( p5| ) for any Pr with t* and preexponent 
factor in (|2^) generally depending on both Pr and Pe, 
this will be the subject of future studies. To conclude, 
the agreement of our result with a variety of data sug- 
gests that Kraichnan model is a proper tool for recovering 
exponents in the theory of turbulent mixing. 
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